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Abstract 

Starting from 2D Euclidean quantum gravity, we show that one recovers 2D Lorentzian 
quantum gravity by removing all baby universes. Using a peeling procedure to 
decompose the discrete, triangulated geometries along a one-dimensional path, we 
explicitly associate with each Euclidean space-time a (generalized) Lorentzian space- 
time. This motivates a map between the parameter spaces of the two theories, un- 
der which their propagators get identified. In two dimensions, Lorentzian quantum 
gravity can therefore be viewed as a "renormalized" version of Euclidean quantum 
gravity. 
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1 Introduction 



The role of baby universes in quantum gravity has been discussed extensively and 
many conjectures have been made concerning their role in providing effective cou- 
pling constants for field theories as well as for gravity itself [|J. Due to the lack 
of a theory of (four-dimensional) quantum gravity, the "integration over baby uni- 
verses" has never been performed in any explicit manner. However, there is a chance 
that this idea can be realized in two dimensions, where quantum gravity is a renor- 
malizable quantum field theory. In this article we will show that starting from 2D 
Euclidean quantum gravity, there is indeed a sense in which one can integrate out 
all baby universes. The resulting "renormalized" quantum theory coincides with the 
theory of so-called Lorentzian 2D quantum gravity. 

The Lorentzian quantum gravity model was introduced in order to have a path- 
integral formulation where the Lorentzian character of the metric is built in at a 
fundamental level. It is also closer to canonical quantization methods, where from 
the outset one works with globally hyperbolic space-times. The two-dimensional 
model can be solved explicitly and is different from Euclidean quantum gravity, by 
which we will mean the theory defined by a path integral over Euclidean geome- 
tries. In the continuum limit, a typical geometry coming from the path integral of 
Lorentzian quantum gravity is two-dimensional E|. This is in contrast with the 
situation in Euclidean quantum gravity, where a typical geometry has an anomalous 
fractal Hausdorff dimension of four 0] . 

We will be using a representation of Euclidean quantum gravity in terms of 
so-called dynamical triangulations. This method provides a reparameterization- 
invariant regularization of the field theory. The two-dimensional geometries will be 
represented by equilateral triangulations (with edges of cut-off length a). The path 
integral is performed by summing over all triangulations of a given topology. The 
continuum limit is obtained as the lattice spacing a —* 0. We will then give an explicit 
prescription for removing baby universes from a given Euclidean triangulation. One 
may view this procedure as analogous to the removal of tadpole graphs in a quantum 
field theory. The latter can usually be viewed as a redefinition of the coupling 
constants and a shift of the fields in the theory. In our case the effect is more 
drastic. We end up with a new, physically inequivalent theory. 

2 The two-loop correlator and the peeling procedure 
2.1 The Euclidean case 

We will investigate the relation between two-dimensional Euclidean and Lorentzian 
quantum gravity by studying the distance-dependent two-loop correlator of the two 
models. The distance between two links I and I' in a given triangulation T is defined 
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as the length of the shortest path connecting I and taken along the links of the 
dual lattice. Furthermore, the distance between a link I and a loop C is defined as 
the minimum distance of the link I to a link belonging to C. With these definitions 
a loop £2 is said to have distance t to a loop C\ if all links in £2 have distance t to 
the loop C\. 

The distance-dependent two-loop correlation function G(li,l2,t) is defined as 

G(h,l 2 ,t)= £ s*CO, (i) 

where the sum goes over all triangulations of cylindrical topology with an exit loop 
£2 of length I2 and an entrance loop £1 of length l\ containing a marked link, the 
distance of £2 from C\ being equal to t. The quantity N(T) is the number of 
triangles in the triangulation T and g is related to the cosmological constant /i by 
g = e-». 

Euclidean and Lorentzian quantum gravity differ by the type of triangulations 
allowed in ([!]). In the Lorentzian case one only sums over triangulations which for 
any t' < t have only one (connected) loop at distance t' from the entrance loop. 
In the usual language of 2D quantum gravity this means that baby-universes are 
forbidden. Furthermore, in the Lorentzian case a loop cannot shrink to length zero. 
Using the distinguished variable t as a time variable, we can divide the links into 
space-like and (future-oriented) time-like, and equip every Lorentzian triangulation 
with a well-defined causal structure By contrast, in the Euclidean case there 
can be more than one loop at a given distance t from the entrance loop and the 
above construction does not work. We note, however, that in || a suggestion of 
introducing a notion of space-like and time-like edges for Euclidean triangulations 
was put forward. This lead to a model equivalent to 2D Euclidean quantum gravity 
coupled to Ising spins. 

It is well-known that a combinatorial identity can be derived for the Euclidean 
version of (|l|) by a so-called peeling procedure §. We will show that also in the 
Lorentzian case a peeling procedure can be defined, which in addition enables us 
to set up a direct correspondence between the discrete Euclidean and Lorentzian 
histories. 

We start by reviewing briefly the peeling formalism in the Euclidean case. In 
order to specify a peeling, one needs a set of moves that systematically remove 
building blocks at a marked link on the boundary of a triangulation and a rule for 
choosing a new marked link on the boundary after each move. We will work in the 
class of so-called unrestricted Euclidean triangulations which apart from triangles 
also contain double links as building blocks. (Unrestricted triangulations constitute 
the class of triangulations appearing in the usual matrix model formulation of 2D 
gravity.) This generalization does not affect any of the universal properties of the 
model. Consider now an unrestricted triangulation contributing to the sum in ([!]). 
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By definition it has an entrance loop with one marked link, which belongs either to 
a triangle or a double link. If it belongs to a triangle we remove the triangle from 
the triangulation. We distinguish three cases, according to whether the triangle 
has one, two or three links in common with the boundary. If it has more than 
one boundary link, additional double links are created during the move, so that the 
length of the original boundary loop is always increased by one (see Fig.l, where 
also the rule for link marking is specified). If the marked link belongs to a double 




Figure 1: The minimal decomposition. In the left column the boundary before the 
decomposition and in the right column the boundary after the decomposition. 

link we simply remove the double link. This will in general lead to a splitting of the 

original loop into two. On each of the resulting loops we place a mark as shown in 

Fig. 1. We can choose either of the two loops as our new entrance loop, associating 

the other one with a baby universe. This gives rise to the factor of two in eqn. (|2|) 

below. (One or both loops may of course be trivial, i.e. consist of only "a marked 

point".) By this procedure we associate with each Euclidean triangulation a unique 

sequence of peeling moves. If the entrance loop has length li, then on average l\ 

peeling steps will move the boundary one time step ahead. Thus one identifies the 

1 d 

deformation 5G(l\,l 2 ,t) associated with a single peeling step with ——G(li,l2,t). 

l\ at 

This deformation fulfils 

—G(l u l 2 ,t) = -l 1 G(l 1 ,l2,t)+ghG(l 1 + l,l 2 ,t) + 2 J2 hW(l)G(h-l-2,l 2 ,t), (2) 
at 1=0 

with the initial condition 

G(h,l 2 ,t = 0) = S hth , (3) 
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where W(l) is the disk amplitude defined by 

W(l)= £ ^ (T) , (4) 
rer(i) 

and the sum is over all (non-restricted) triangulations with one connected boundary 
component of length I. Let us introduce the generating functionals 

oo 1 1 

G(z,y,t) = £ G (! lM) ( 5 ) 

h,h=o z y 

oo 

w ^ z ) = E— ^(0- (6) 

1=0 z 

Here z and y can be understood as boundary cosmological constants. Then the 
differential equation (|2|) reads 

G(z,y,t) = — (h(z)G(z,y,t)), (7) 
at oz 



with initial condition 



where 



G(z,y,t = 0) = — — !— -, (8) 
zy zy - 1 



= z - gz z - 2W(z). (9) 

It is well known that in analogy with the procedure described above, it is possible 
to derive a combinatorial identity for W(z) by removing a triangle or a double 
link from a triangulation contributing to the sum in (||), see for instance Q. The 
combinatorial identity reads 

(z - gz 2 )W(z) - 1 + g{ Wl {g) + z) = (W(z)f , (10) 

where 

M9)= E 9 N{T) - (11) 
TeT(i) 

Here the terms proportional to g emerge when a triangle is removed whereas the 
term (W(z)) 2 emerges when a double link is removed and the surface splits into two. 
We note that if we set g = we get 

W(z) = \{z- ^~4) , (12) 

which is the well-known generating function for rooted branched polymers. These 
appear exactly because of our inclusion of double links. For g ^ 0, VF(z) is the 
solution of a quadratic equation and the unknown constant w\ (g) is determined by 
noticing that the analyticity structure of W(z) should not change discontinuously 
at g = 0. The solution for W(z) reads 

W(z) = l -(z-g Z i + g( z - c (g)) yj (z - c + (g) ) (z - c_ (g) )) , (13) 
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where c—(g) < c+(g) < c(g) and where all three quantities can be expressed explicitly 
in terms of g. At a certain value g c of g the radius of convergence of the series @ 
is reached and a continuum limit can be defined. Setting 

g c - g = a 2 A, (14) 

where a? is a scaling parameter with the dimension of area and A is the continuum 
cosmological constant we have for g close to g c 

c(g) = z c (l + ±ay/A) +0(a 2 ), (15) 

c+(g) = z c (l-aVX) +0(a 2 ), (16) 
c_( 5 ) = c_+0(a 2 ). (17) 

It is now natural to introduce a continuum boundary cosmological constant Z by 

z = z c (l + aZ). (18) 

We thus have that h(z) exactly consists of the scaling part of —2W(z). Introducing 
a continuum time variable T by 

T = z]J\z c - c^fl 2 g c a x l\ (19) 

one obtains a differential equation for the continuum version G(Z, Y, T) of the two- 
loop correlator, namely, 

±G(Z,Y,T) = --^ ^(Z - ±VX)y/z + VA) G(Z,Y,T)j , (20) 
with the initial condition 

G(Z,Y,T) = T L p . (21) 

Note that the continuum cosmological constants A and Z and the time T are only 
defined up to positive multiplicative constants. We have used this freedom to bring 
the propagator equation into the simple form (]2(i|). 

2.2 The Lorentzian case 

One may invert the peeling procedure introduced in the previous section to obtain 
a Euclidean triangulation from a sequence of peeling moves. Its geometry is given 
in the form of a bulk cylindrical geometry connecting the entrance and exit loops, 
"decorated" with baby universe out-growths. This is the case because some of the 
peeling steps remove just a single triangle (which we now associate with the bulk 
cylinder), while the removal of a double link in general amounts to the deletion of 
an entire baby universe. 

As mentioned earlier, the difference between Euclidean and Lorentzian gravity 
can be traced to the presence or absence of baby universes. There is already an 



5 



explicit construction of Euclidean from Lorentzian histories through the addition of 
baby universes j^]. Our current aim is to go the other way. The peeling procedure 
suggests an obvious map from Euclidean to Lorentzian discrete geometries. Given 
a representation of a Euclidean geometry as a sequence of peeling moves, consider a 
step where a double link-cum-baby universe is removed, see Fig. 2. Decompose the 



Figure 2: A representation of a Euclidean triangulation as a sequence of peeling 
moves along a "spiral-like" one-dimensional curve. Black rectangles are double links 
and grey circles baby universes. 

baby universe into a branched polymer attached to the marked link and a smaller 
baby universe attached to the branched polymer in such a way that the branched 
polymer consists of the largest possible number of links. Next, replace the step 
where the double link-cum-baby universe is removed by a step where only the double 
link and the branched polymer is taken way. Do the same for all similar steps in 
the sequence of peeling moves. The resulting sequence of peeling moves obviously 
represents a geometry without baby universes. It is clear that the class of geometries 
one generates this way is larger than the class of Lorentzian triangulations introduced 
in Q . Firstly, it contains all of the regular Lorentzian triangulations of Q , as can be 
checked explicitly. Secondly, the boundaries may now contain double links. We will 
refer to this generalized class of geometries as "Lorentzian geometries". This will 
be justified below by showing that they lead to a theory in the same universality 
class as the original Lorentzian model and the extensions discussed in Q. The 
differential equation for the two-loop correlator of Lorentzian geometries is obtained 
by replacing 2W(z) in (||) by Wbp(z), and thus reads 

d „ ~, d 



d -M(~z, y, t) = — (h(z)G(z, y, t)) , (22) 



with initial condition 



where 



G(2,jf,t = 0) = i— (23) 
zyl- zy 



h{~z) = z-g~z 2 - W B p{z). (24) 



We have introduced tildes on all the Lorentzian variables in order to distinguish 
them from the Euclidean ones Let us note here that it is easy to generalize this 



6 



model of Lorentzian quantum gravity to models with arbitrary polygons as building 
blocks. Using the newly developed technology of || it may be possible to solve all of 
these models exactly. Here we only consider the simplest model. We wish to locate 
a critical point (jg c , z c ) at which we can define a non-trivial continuum limit by the 
scaling relations 

g = g c (l - cia 2 A), z = z c (l + c 2 aZ) , (25) 

where c\ and c 2 are numerical constants. The only way to obtain a non-trivial 
differential equation in the continuum limit is by demanding that h(z) ~ 0(a 2 ). 
This determines the critical values of z and g to be 

z c = 7f! g c = — . (26) 

For later convenience we set c 2 = and c\ = -^=. We also introduce a continuum 
version T of t by 

t = -f. (27) 
a 

Note that the relative dimensions of the time variable and the area are different for 
Euclidean and Lorentzian quantum gravity. Then we derive the following evolution 
equation for the continuum two-loop correlator G(Z, Y, T) 

■^G(Z,Y,f) = --t ((Z 2 - A) G(Z,Y,fj) , (28) 
with initial condition 

G(Z,Y,f) = ^J-^. (29) 
Z + Y 

This differential equation coincides with the one derived in using a completely 
different strategy. 

Let us now turn to a discussion of the disk amplitude W(l) of Lorentzian gravity. 
At the discrete level it may be defined as the sum J2?Li h = 0, t). This reflects 
the fact that a Lorentzian geometry with a regular initial spatial geometry at t = 
can only terminate at a finite t if we allow for a singularity where space contracts 
to a point. 

We can derive a combinatorial identity for the generating functional W{z) in the 
same way as for the Euclidean disk amplitude W(z). The combinatorial identity 
reads 

(~z - gz 2 )W(z) - 1 + giMg) + *) = W BP (z)W(z). (30) 

As in the case of the two-loop correlator, the only modification of the equation 
concerns the term associated with the removal of a double link. If, in the Euclidean 
case, the marked link is part of a double link, we remove the double link and are 
left with two disjoint disks. If, in the Lorentzian case, the marked link is part of a 
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double link, we remove the double link and are left with one disk and one branched 
polymer. In this case, by setting g = 0, we get 

W{z) = W BP (~z). (31) 

This underlines the consistency of our construction. For g 7^ on the other hand 
we obtain 

t(^ 2( '-" (32) 
z - 2gz 2 + Vz 2 - 4 

Note that the denominator of (55) is equal to 2h(z). For g < \ this function has one 
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zero zi(g), which behaves as i + O(g) when g is close to zero. Hence a pole appears 



in the transition from (31) to ( p2|) and we must choose the unknown constant w\{g) 
such that this pole is cancelled. At g = 4 the function h(z) acquires an additional 
zero, Z2(g), and the critical point g c = corresponds to the situation where the 
two zeros z\(g) and £2(5) coincide. In the vicinity of the critical point we have the 
expansions 

4 16 

h(g) = — + - J -(g c -g)W, (33) 

M9) = -^-^(~g c -g) 1/2 . (34) 

Using (33), fl34|) and p5|), we derive the continuum disk amplitude 

W(Z) = lim aW(z) = -L^-^. (35) 
a^o 3 5 / 4 z + V A 

This expression coincides with the one found in P] by completely different means. 



3 Renormalizing Euclidean gravity to Lorentzian grav- 
ity 

We will now show that it is possible to "renormalize" the Euclidean coupling con- 
stants and the Euclidean time such that the equation governing the time evolution 
of the Euclidean two-loop function is mapped to the corresponding equation for the 
Lorentzian two-loop function. 

We start from the discrete version of the differential equation for the time devel- 
opment of the Euclidean universes, (0), and make an ansatz 

z = z{z,g), *=T*> 9 = 9(9), (36) 

where b may depend on g and g, but not on z or z. In principle it should be possible 
to determine the required map between the coupling constants of the two theories 
directly from the peeling procedure which maps Euclidean to Lorentzian histories. 
However, it turns out to be simpler to use an indirect argument. // there exists a 
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transformation of the form ( pq ) which maps the Euclidean differential equation into 
the Lorentzian differential equation, then clearly it must include a wave- function 
renormalization of the two-loop correlator 



and furthermore 



or, equivalently, 



G(z,y,t) oc ( ■£) G(z(z),y(y),bt), 



bh(z)—=h(z) 



dz r dz 

b I ^— + k, 



(37) 

(38) 
(39) 



h(z) J h (z) 
where k is an integration constant. 

Using the explicit expressions @ and ( p4[ ) for h(z) and h(z), we can carry out 
the integrations in (^) and find 

1 J 25(g) + e(g)(z-c(g)) + 2^5(g)^(z-c + (g))(z-c.(g)) \ 

VW) I J 1 j 

log(z - zi(#)) + — 7 — rr-log(z - z 2 (s)) 



2? 



2~z x (g) [zM- | 



25 2 (ff) (SzCfir) - | 



zi(g) 2 - 4\/5 2 - 4 + zz\(g) - 4 



+ ... 



+ fc, 



Here . . . means terms similar to the previous one with z\(g) replaced by 2:2(5), 23(g) 
and Zi(g) where z\(g), ^(g), z%(g) and z^(g) are the four roots of the polynomial 



r 4 1 ~3 1 1 

a g 2 

page. Furthermore 5 and e are defined by 



: ' - l z + jj with 21(5) and 22(5) still playing the roles described on the previous 



5(g) = (c(g)-c + (g))(c(g)-c^(g)), 

e (s0 = 2c (s0 - c +(s0 - c -(fiO- 



For any value of k, g, g and b, the relation ( |4CJ| ) defines a mapping 5 = i(z). However, 
not all values of b and k correspond to situations we can encounter in the peeling 
procedure. Firstly, the two theories have to approach criticality simultaneously, 
since macroscopic surfaces and boundaries are created at the critical points. Sec- 
ondly, we insist on a canonical scaling dimension of the cosmological and boundary 
cosmological constants. This amounts to saying that the scaling of the continuum 
cosmological constants is governed by the same power of the lattice cut-off a. Since 
we have assumed that g is only a function of g, and not of z and z, this implies 
that A oc A, with a proportionality coefficient independent of the cut-off. In the 
scaling limit all dependence on z^(g) and z^(g) disappears and using the scaling 
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relations (|i^~(|lq), (pq), (|33;) and (|34j), we see that in order for the mapping to 
persist in the continuum limit, b and k must scale according to 

b = a 1 ' 2 A 1 / 4 /?, fe = o- 1 / 2 A- 1 \ (41) 



for two dimensionless constants /3 and k. To fix the relevant power of A in (41), we 
have used the fact that b and k cannot depend on z or z. Note that b and k behave 
in a singular way when we approach the continuum limit. This is an unavoidable 
consequence of interpolating between two universality classes of 2D geometries with 
different scaling behaviour (recall that in the previous section we showed that t scales 
as t ~ a -1//2 , whereas t ~ a -1 ). 

It is now straightforward to compute the continuum version of (fi"0|), 



-^log XI V , = ^log + — , (42) 




3 

where we have absorbed some additional numerical constants into j3 and k. We will 
argue that in order for this map to be associated with the cutting procedure we must 
choose k = 0. Taking the limit Z —* oo implies that the boundary length L —* 0. 
Now, a Lorentzian boundary is shorter than a Euclidean one, since it was obtained 
from the latter by cutting away baby universes which at a given time contributed 
to the boundary length. This is the physical motivation for requiring that Z — > oo 
should imply Z — * oo, i.e. k = 0. Moreover, a rescaling of both Z and y/X with a 
common factor allows us to choose 

VA = ivA> (43) 
in such a way that the map (fS) simplifies to 



+ (44) 



V3 AV ' 

and the relation between the continuum times T and T becomes 

T 



/3AV4- 



(45) 



One can of course verify directly that the (non-analytic) change of variables (03) and 
( |45| ) supplemented by the wave function renormalization ( |37| ) transforms eq. ( p0[ ) 
for Euclidean quantum gravity into eq. (|28| ) for Lorentzian gravity^]. To map the 
actual two- loop correlator of Euclidean quantum gravity Go(Z, Y, T) (i.e. the solution 
to (|2"C)|) with initial condition (|2l"D) onto the actual two-loop correlator of Lorentzian 



Note that some additional constants have been absorbed into /3 in going from © to @ 
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quantum gravity Gq(X,Y,T) (i.e. the solution to ( p8[ ) with initial condition (p9|)) 
an additional dressing factor is needed More precisely one has 

G (Z,Y,f) = K(Z(T,Z),Y) G (Z,Y,T), (46) 

where 



K(S(T,Z),Y)= - mn+Z, / Z + ^ , (47) 

Z(Z(T,Z))+Y(Y) Jz(T,Z) + SK 



and where Z(T,Z) is the solution to characteristic equation for (|20|): 

=-(Z-1VA) V / ^ + v / A, Z(T = 0,Z) = Z. (48) 

Summing up, we have shown that "integrating over baby" universes can be described 
as a renormalization of the cosmological constants and the time variable combined 
with a dressing of the two-loop correlator. 



4 Discussion 



We have shown how one can analyze the geometric structure of Lorentzian and Eu- 
clidean triangulations by representing them as a sequence of peeling moves defining 
a one-dimensional 'spiral-like' path. Along this path, one meets in both cases ei- 
ther triangles or double links. The difference between the Euclidean and Lorentzian 
cases is that in the former a double link may have a branched polymer plus an entire 
baby universe attached to it whereas in the latter a double link can have at most 
a branched polymer attached to it. At the level of sequences of peeling moves, this 
allows us to define a many-to-one map from Euclidean to Lorentzian geometries. 
This is reminiscent of the situation one encounters in analyzing the fractal structure 
of random walks versus branched polymers. There, moving along a path connecting 
two vertices, in the random walk case one may have only a single link emerging at 
each vertex. For the branched polymer case, there can be an entire polymer attached 
at each vertex, see e.g. 0. In this context it may be of interest to point out that 
pure Lorentzian quantum gravity has been proven to be equivalent to a model of 
random walks []|]. 

This analogy is further strengthened by the nature of the continuum renormal- 
ization of the coupling constants presented in section ||. Apart from trivial rescal- 
ings, only the boundary cosmological constant is subject to a renormalization. This 
reflects the fact that the peeling decomposition proceeds along a one-dimensional 
path. Moreover, the functional form of the relation between the Lorentzian and 
Euclidean boundary cosmological constants is very similar to that between the cos- 
mological constants of random walks and branched polymers, since in both cases 
the cosmological constant of the simpler model has a square-root dependence on the 
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cosmological constant of the more extended model. Furthermore, in both cases the 
fractal dimension of the decorated system is twice that of the simpler system. 
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